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Abstract. We introduce some multiple integrals that are expected to have 
the same singularities as the singularities of the n-particle contributions \^ 
to the susceptibility of the square lattice Ising model. We find the Fuchsian 
linear differential equation satisfied by these multiple integrals for n = 1, 2, 3,4 
and only modulo some primes for n = 5 and 6, thus providing a large set of 
(possible) new singularities of the x ■ We discuss the singularity structure 
for these multiple integrals by solving the Landau conditions. We find that the 
singularities of the associated ODEs identify (up to n = 6) with the leading 
pinch Landau singularities. The second remarkable obtained feature is that the 
singularities of the ODEs associated with the multiple integrals reduce to the 
singularities of the ODEs associated with a finite number of one dimensional 
integrals. Among the singularities found, we underline the fact that the quadratic 
polynomial condition 1 + 3ui + iw 2 = 0, that occurs in the linear differential 
equation of x' 3 '' j actually corresponds to a remarkable property of selected elliptic 
curves, namely the occurrence of complex multiplication. The interpretation of 
complex multiplication for elliptic curves as complex fixed points of the selected 
generators of the renormalization group, namely isogenies of elliptic curves, is 
sketched. Most of the other singularities occurring in our multiple integrals are 
not related to complex multiplication situations, suggesting an interpretation in 
terms of (motivic) mathematical structures beyond the theory of elliptic curves. 



PACS: 05.50. +q, 05.10.-a, 02.30.Hq, 02.30. Gp, 02.40.Xx 

AMS Classification scheme numbers: 34M55, 47E05, 81Qxx, 32G34, 34Lxx, 
34Mxx, 14Kxx 

Key-words: Susceptibility of the Ising model, singular behaviour, Fuchsian 
linear differential equations, apparent singularities, Landau singularities, pinch 
singularities, modular forms, Landen transformation, isogenies of elliptic curves, 
complex multiplication, Heegner numbers, moduli space of curves, pointed curves. 

1. Introduction 

The susceptibility x °f the square lattice Ising model has been shown by Wu, McCoy, 
Tracy and Barouch [I] to be expressible as an infinite sum of holomorphic functions, 
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given as multiple integrals, denoted X ; that is kT ■ \ — ■ B. Nickel 

found [2j [3] that each of these x^' s is actually singular on a set of points located 
on the unit circle \s\ = |sinh(2i^)| = 1, where K = J/kT is the usual Ising model 
temperature variable. 

These singularities are located at solution points of the following equations: 



(• + ;) 



IN ft 1 m 1 

u H r + u m H 



u 2n+1 = 1, -n < to, k < n (1) 

From now on, we will call these singularities of the "Nickelian type", or simply 
"Nickelian singularities". The accumulation of this infinite set of singularities of the 
higher-particle components of x( s ) on the unit circle \s\ — 1, leads, in the absence 
of mutual cancellation, to some consequences regarding the non holonomic (non D- 
finite) character of the susceptibility, possibly building a natural boundary for the 
total x( s )- However, it should be noted that new singularities that are not of the 
"Nickelian type" were discovered as singularities of the Fuchsian linear differential 
equation associated [H[H[6] with x^ and as singularities of x^ itself [7] but seen as 
a function of s. They correspond to the quadratic polynomial 1 + 3 w + 4u> 2 where 
2w = s/(l + s 2 ). In contrast with this situation, the Fuchsian linear differential 
equation, associated [5] with x' 4 \ does not provide any new singularities. 

Some remarkable Russian-doll structure as well as direct sum decompositions were 
found for the corresponding linear differential operators for x*' 3 ' and x*' 4 '- In order 
to understand the "true nature" of the susceptibility of the square lattice Ising model, 
it is of fundamental importance to have a better understanding of the singularity 
structure of the n-particle contributions x^ , and also of the mathematical structures 
associated with these x > namely the infinite set of (probably Fuchsian) linear 
differential equations associated with this infinite set of holonomic functions. Finding 
more Fuchsian linear differential equations having the x^'s as solutions, beyond those 
already found [H [8] for x' 3 ' and x' 4 ' ; probably requires the performance of a large 
set of analytical, mathematical and computer programming "tours-de-force". 

As an alternative, and in order to bypass this "temporary" obstruction, we have 
developed, in parallel, a new strategy. 

We have introduced [7] some single (or multiple) "model" integrals as an "ersatz" 
for the x^' s as far as the locus of the singularities is concerned. The x^' s are 
defined by (n — l)-dimensional integrals OmUOj (omitting the prefactoiQ) 



*<») = » 
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where 



and 



j — 1 l<i<j<n •* 



Ri n) = 1 + nti *i (4) 



i - n: =1 

t The prefactor reads (1 - s 4 ) 1 / 4 /* for T > T c and (1 - s" 4 ) 1 / 4 for T < T< 
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with 

2w 



1 - + ^(1 - 2w cos(0,)) 2 - Aw 2 



(5) 



ft = y X> = 0. (6) 

y (1 - 2wcos(0i)) - 4w 2 3=1 

The two families of integrals we have considered in [7] are very rough 
approximations of the integrals |(2j). For the first familjjj], we considered the n-fold 
integrals corresponding to the product of (the squar^l of the) UiS, integrated over the 
whole domain of integration of the 0, (thus getting rid of the factors G* 1 ' and R^). 
Here, we found a subset of singularities occurring in the \ as well as the quadratic 
polynomial condition 1 + 3w + Aw 2 = 0. 

For the second family, we discarded the factor G^ and the product of j/j's, 
and we restricted the domain of integration to the principal diagonal of the angles 
0i (01 = 02 = ••• = 0>i-i )• These simple integrals (over a single variable), were 
denoted [7] 

#w = -i +i pa* i ( 7 ) 

n! n\ J 2tt 1 -x n - l ((j>) -x((n- 1)0) W 

where £c(0) is given by ©. 

Remarkably these very simple integrals both reproduce all the singularities, 
discussed by Nickel [21 E], as well as the quadratic roots of 1 + 3w + Aw 2 = 
found [IIH] for the linear ODE of x' 3 -*. One should however note that, in contrast 
with the x*") , no Russian-doll or direct sum decomposition structure is found for the 

(n) 

linear differential operators corresponding to these & D . 

Another approach has been introduced as a simplification of the susceptibility 
of the Ising model by considering a magnetic field restricted to one diagonal of the 
square lattice [11] , For this "diagonal susceptibility" model [IT], we benefited from 
the form factor decomposition of the diagonal two-point correlations C(N,N), that 
has been recently presented [12], and subsequently proved by Lyberg and McCoy [13j. 

(n) 

The corresponding n-fold integrals x d were found to exhibit remarkable direct sum 
structures inherited from the direct sum structures of the form factor |11[ [T2] . The 
linear differential operators of the form factor |12j being closely linked to the second 
order differential operator Le (resp. Lk) of the complete elliptic integrals E (resp. 
K), this "diagonal susceptibility" model pj] is closely linked to the elliptic curves of 
the two-dimensional Ising model. By way of contrast, we note that the singularities of 
the linear ODE's for these n-fold integrals [11] are quite elementary (consisting 
of only n-th roots of unity) in comparison with the singularities we encounter for the 
integrals on a single variable |(7|). 

These two approaches corresponding to two different sets of n-fold integrals of 
the Ising class [14] are complementary: ([7]) is more dedicated to reproduce the non- 
trivial head polynomials encoding the location of the singularities of the j but fails 

t Denoted Y^(w) in [7]. 

| Surprisingly the integrand with (OJ-i%) 2 yi e lds second order linear differential equations 0, 
and consequently, we have been able to totally decipher the corresponding singularity structure. By 
way of contrast the integrand with the simple product (FJ fjj) yields linear differential equations 
of higher order, but with identical singularities 0- 
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to reproduce some remarkable (Russian-doll, direct sum decomposition) algebraico- 
differential structures of the corresponding linear differential operators, while the other 
one [TT] preserves these non-trivial structures of the corresponding linear differential 
operators but provides a poorer representation of the location of the singularities (n-th 
roots of unity) . 

In this paper, we return to the integrals ^ where, this time, the natural next 
step is to consider the following family of n-fold integrals 



which amounts to getting rid of the (fermionic) factor (G^) 2 in the n-fold integral 
(f5|). This family is as close as possible to (J2j) , for which we know that finding the 
corresponding linear differential ODE's is a huge task. The idea here is that the 
methods and techniques we have developed [4j [5] for series expansions calculations 
of x and x^\ seem to indicate that the quite involved fermionic term (G^) 2 in 
the integrand of ^ should not impact greatly on the location of singularities of these 
n-fold integrals (J2]) ■ This is the best simplification of the integrand of <[2j) for which we 
can expect to retain much exact information about the location of the singularities of 
the original Ising problem. However, we certainly do not expect to recover from the 
n-fold integrals ([8]) the local singular behavior (exponents, amplitudes of singularities, 
etc ...). Getting rid of the (fermionic) factor (G^) 2 are we moving away from the 
elliptic curves of the two-dimensional Ising model ? Could it be possible that we lose 
the strong (Russian-doll, direct sum decomposition) algebraico-differential structures 
of the corresponding linear differential operators inherited from the second order 
differential operator Le (resp. Lk) of the complete elliptic integrals E (resp. K) 
but keep some characterization of elliptic curves through more "primitive" (universal) 
features of these n-fold integral like the location of their singularities ? 

In the sequel, we give the expressions of , and the Fuchsian linear 

(n) 

differential equations for & H for n = 3 and n = 4. For n = 5, 6, the computation 
(linear ODE search of a series) becomes much harder. Consequently we use a modulo 
prime method to obtain the form of the corresponding linear ODE with totally explicit 
singularity structure. These results provide a large set of "candidate singularities" for 
the ■ From the resolution of the Landau conditions [7] for (JSj) , we show that the 
singularities of (the linear ODEs of) these multiple integrals actually reduce to the 
concatanation of the singularities of (the linear ODEs of) a set of one-dimensional 
integrals. We discuss the mathematical, as well as physical, interpretation of these 
new singularities. In particular we will see that they correspond to pinched Landau- 
like singularities as previously noticed by Nickel |15j . Among all these polynomial 
singularities, the quadratic numbers l+3w+Aw 2 — are highly selected. We will 
show that these selected quadratic numbers are related to complex multiplication for 
the elliptic curves parameterizing the square Ising model. 

The paper is organized as follows. Section ^ presents the multidimensional 

(n) 

integrals Q H and the singularities of the corresponding linear ODE for n = 3, ■ • ■ , 6 , 
that we compare with the singularities obtained from the Landau conditions. We 
show that the set of singularities associated with the ODEs of the multiple integrals 

(n) 

$ H reduce to the singularities of the ODEs associated with a finite number of one- 
dimensional integrals. Section ((3j) deals with the complex multiplication for the elliptic 
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curves related to the singularities given by the zeros of the quadratic polynomial 
1 + 3w + 4 w 2 = 0. Our conclusions are given in section ([4]). 

2. The singularities of the linear ODE for 

For the first two values of n, one obtains 

and 



= 2 ' l-i 6w 2 ■ ^i(l/2, -1/2; 1; W). (10) 

(n) 

For n > 3, the series coefficients of the multiple integrals § H are obtained by 
expanding in the variables Xi and performing the integration (see |Appendix A[ ) . One 
obtains 



1 1 



X 



= 1 ' EE ( 2 - **.o) • ( 2 - <W) ' ™ n(k+p) ■ « n (M (11) 

fe=0p=0 

where a(k,p) is a 4F3 hypergeometric series dependent on w. 

The advantage of using these simplified integrals JHJ instead of the original ones 
|2j is twofold. 

Using (flT|) the series generation is straightforward compared to the complexity 
related to the x^ . As an illustration note that on a desk computer, $jj are generated 
up to w 200 in less than 10 seconds CPU time for all values of n, while the simplest 
case of the X > namely x*- 3 -*, tooks three minutes to generate the series up to w 200 . 
This difference between the $^ and x^ increases rapidly with increasing n and 

(n) 

increasing number of generated terms. We note that for the $> H quantities and for a 
fixed order, the CPU time is decreasing with increasing n. For x^ the opposite is 
the case. 

The second point is that, for a given n, the linear ODE can be found with less 
terms in the series compared to the linear ODE for the x ■ Indeed for x^ 3 \ 360 

(3) 

terms were needed while 150 terms were enough for &x . The 

same feature holds for 

X (4) and $^ (185 terms for x (4) and 56 termsfij for $^ } ). 

With the fully integrated sum (ffTjl . a sufficient number of terms is generated 
to obtain the linear differential equations. We succeeded in obtaining the linear 
differential equations, respectively of minimal order five and six, corresponding to 

and . These linear ODE's are given in | Appendix B| 



For $^ (n > 5), the calculations, in order to get the linear ODEs become really 
hugdj. For this reason, we introduce a modular strategy which amounts to generating 
large series modulo a prime and then deducing the ODE modulo that prime. Note that 
the ODE of minimal order is not necessarily the simplest one as far as the required 
number of terms in the series expansion to find the linear ODE is concerned. We have 
already encountered such a situation (HUH]. For (resp. $^ '), the linear ODE 

tt This can be seen from the series expansion {TTJ. Denoting -Ro the fixed order, one has n- (p + fc) < 
-Ro, while the CPU time for the series generation of a n (k, p) is not strongly dependent on n. 
1f From now on, for even n, the number of terms stands for the number of terms in the variable 
x = w 2 . 

X Except the generation of large series which remains reasonable. 
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of minimal order is of order 17 (resp. 27) and needs 8471 (resp. 9272) terms in the 
series expansion to be found. 

Actually, for $^ (resp. $#''), we have found the corresponding linear ODEs of 
order 28 (resp. 42) with only 2208 (resp. 1838) terms from which we have deduced 
the minimal ones. 

The form of these two minimal order linear ODEs obtained modulo primes is 
sketched in |Appendix B| In particular, the singularities (given by the roots of the head 
polynomial in front of the highest order derivative) , are given with the corresponding 
multiplicity in |Appcndix B| Some details about the ODE search are also given in 
|Appendix B 

We have also obtained very long series (20000 coefficients) modulo primes for 
® H , but, unfortunately, this has not been sufficient to identify the linear ODE (mod. 
prime) up to order 100. 

The singularities of the linear ODE for the first are respectively zeros of the 
following polynomials (besides w = oo): 

n = 3, to - (l - 16w 2 ) (1 - to) (1 + 2 to) (l + Zw + Aw 2 ) , 

n = 4, to ■ (l - 16 w 2 ) (l - Aw 2 ) , 

n = 5, to • (l - 16 w 2 ) (1 - w 2 ) (1 + 2 w) (l + 3 to + 4 w 2 ) 
(l~3w + w 2 ) (l + 2w-4w 2 ) (l + 4w + 8w 2 ) 
(l - 7w + 5w 2 - Aw 3 ) (l-w-3w 2 + 4io 3 ) 
(l + 8w + 20w 2 + 15w 3 + 4to 4 ), (12) 

n = 6, to - (l- 16w 2 ) (1 -Aw 2 ) (l - w 2 ) (l - 25w 2 ) 
(1 - 9 w 2 ) (1 + 3w + Aw 2 ) (l-3u; + Aw 2 ) 
(l-10w 2 + 29w 4 ). (13) 

For n — 7 and n = 8, besides modulo primes series calculations described above, 
we also generated very large series from which we obtained in floating point form, the 
polynomials given in |Appendix C| (using generalised differential Pade methods) . 

If we compare the singularities for $^ to those obtained with the "Diagonal 
modeU presented in [7], i.e. one sees that the singularities of the linear ODE 

(n) 

for the "Diagonal model" are identical to those of the linear ODE of the $> H for 
n = 3, 4 (and are a proper subset to those of $^ for n = 5, 6). The additional 
singularities for n = 5,6 are zeros of the polynomials: 

n = 5, (l + 3w + Aw 2 ) (l + Aw + 8w 2 ) 

x (l - 7w + 5w 2 - Aw 3 ) , 
n = 6, (l +3w + Aw 2 ) (l - 3to + Aw 2 ) (l - 25 w 2 ) . 

For n = 7, the zeros of the following polynomials (among others) are singularities 
which are not of Nickel's type |T|) and do not occur for 

1 + 8w + I5w 2 - 2lw 3 - 60 w A + 16 w 5 + 96 w 6 + 6Aw 7 , 
1 - Aw - IQw 2 - 48w 3 + 32 w 4 - 128 w 5 . 
tt Not to be confused with the "diagonal susceptibility" and the corresponding n-fold integrals 



Singularities of n-fold integrals 



7 



The linear ODEs of the multiple integrals & H thus display additional 
singularities for n = 5, 6 and n = 7 (n = 8 see below) compared to the linear 
ODE of the single integrals 3>jj . 

We found it remarkable that the linear ODEs for the integrals 4>£^ display all the 
Nickelian singularities, as well as the new quadratic numbers 1 + 3w + Aw 2 = found 
for ■ It is thus interesting to see how the singularities for <£>^ are included in 
the singularities for and whether the new (with respect to singularities 
can be given by one-dimensional integrals similar to & D . Let us mention that 
the singularities of the linear ODE for (respectively ) are remarkably also 

singularities of the linear ODE for $^ (respectively ^^)- ^ n t° e following, we will 
show how this comes about and how it generalizes. For this, we solve in the sequel 
the Landau conditions for the n-fold integrals ([5]). 

(n) 

2.1. Landau conditions for the $> H 

We remind the reader that the Landau conditions [7] are necessary conditions for 
singularities to be the singularities of the integral representation itself. In a previous 
paper [7], we have shown for particular integral representations belonging to the 
Ising class integrals [14], that in fact the solutions of Landau conditions identify for 
specifi(0 configurations (see below) with the singularities of the ODE associated with 
the quantity under consideration. 

The Landau conditions [7] amount to carrying out algebraic calculations [7] on 
the integrand ((Sj) to get singularities of these n-fold integrals or even, as we will see 
in the sequel, singularities of the corresponding linear ODE [7]. 

In the sequel we use the following integral representation (HE]: 




(14) 



Defining 

D{<pj,tpj) = 1 — 2w ■ (cos(4>j) + cos(^i)) , 



(15) 



the integral (see its expansion IjA.ip in | Appendix A[ ), becomes 




(16) 



where the Dirac delta's are introduced to take care of the conditions 



n 





mod. 2 7r 



(17) 



on both the angles 4>j and the auxiliary angles ipj. 



% In that respect one must recall the notion of leading singularities in contrast with the subleading 
singularities (see page 54 in [16].) 
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The Landau conditions [HU [17] can easily be written [7] : 

atj ■ D ((j> 3 ,ipj) = 0, j = 1, • • • , n, (18) 

Pj'4>j=0, 7r^=°: j= 1, n-1, (19) 

ckj • sin(#j) - a„ • sin(0„) + = 0, j = 1, • • • n - 1, (20) 

ay • sin(-0j) - a„ • sin(^„) + 7,- = 0, j = 1, • • • , n ~ 1 (21) 

together with (fT7|) . The Landau singularities are obtained by solving these equation^ 
in all the unknowns, where the parameters ctj, (3j, jj should not be all equal to zero. 
In this paper, our aim is not to find all the solutions of the above equations but to 

(n) 

show that the singularities of the linear ODE for the & H are solutions of the Landau 
conditions. Furthermore, in working out various Ising class integrals [I4j and the two 
models of [7] (see | Appendix D |, we remarked that the singularities of the linear ODE 



are, in fact, included in a particular "configuration". What we mean by "configuration" 
is the set of values (equal to zero or not) of the parameters ctj, (3j, jj, 
The "configuration" we consider 

a 3 + 0, ft = 7j = 0, (22) 

corresponds to pinch singularities on the manifolds D (ipj, ipj) = 0. One may also be 
convinced to take {3j = jj = 0, since the integrand is periodicHj] in <j)j and ipj- 

Let us stress that the the configuration considered where all the Lagrange 
multipliers of the singularity manifolds D(4>, tp) are different from zero (ctj ^ 0, for 
any j) leads to the so-called leading Landau singularities following the terminology of 
page 54 of [16J. 

The Landau conditions become: 

1 - 2 w ■ (cos(0j) + cosiipj)) =0, j = 1, ■ ■■ , n, (23) 

ctj sin(^j) - a n sin(0 n ) =0, j = 1, • • • , n — 1, (24) 

ctj sin(^)j) - a n sm(ip n ) =0, j = 1, • • • , n - 1. (25) 

and: 



E 



= 0, ■ = mod - 2 n ( 26 ) 

j=i i=i 
The Landau singularities are solutions of these conditions (see |Appendix E| for 
details). Note that the first three conditions (|23"|) . (f24|) . ([2511 are invariant by the 
transformation: 

W ► — W, ► Ipj + TT, Ipj > Ipj + TT. (27) 

but the Landau conditions |23|) . |24|) . l[25|) together with ([26)1 are invariant by 
transformation |27j) if and only if n is even. This distinction between even and 
odd integer n (corresponding to the symmetry breaking of w «-> —w) is reminiscent 
of the distinction between even and odd integer n for the associated with the 
distinction between low and high temperature regimes. 

The Landau conditions yield two families of singularities expressed in terms of 
Chebyshev polynomials of the first and second kind. The first family reads: 

T 2pi (1/2W+1) = T n - 2pi -2 P2 (l/2w - 1) , (28) 

< pi < [n/2], < p 2 < [n/2] - Pl 

jt Note that conditions l|19|l . /3j ■ <l>j = 0, jj ■ = 0, j = 1, ■ • • , n — 1 have to be considered in the 
general Landau conditions. They do not occur if one restricts oneself to pinch singularities, 
ff B. Nickel, private communication. 
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The second family is given by the elimination of z from: 

4u> — z 



(Quo — z \ 







\2w Z ) 


~Tn 2 ( 





1 Aw — z 



2w 1 — iw z 



= 0, 



- r,„. , l — — : ) -r,,, t ( z ) = o 



4w z 



with 



?;i = pi, n 2 = n - pi -2p 2 , (30) 

< Pi < n, < p 2 < [{n- Pl )/2]. (31) 

One recognizes in the first set of equations (|28| . a generalization of the 
singularities given by Nickel [15] for the pinch singularities coming from the product 
of the yj's, and also derived for our multiple integral denoted in [7]. These have 
been written as [7} [15]: 

T k (l/2w + l) = T n _ fc (l/2tu - 1) (32) 

Note that, comparatively to (|28)l . the integer fc should be eve 

The second set of equations ([29]) is a generalization of the singularities we derived 
for $^ in [7]. In both formulae, one notes the occurrence of a second varying integer 
P2, leading to a better understanding of the singularities of these integrals. Indeed with 
P2 running, the linear ODE for & H will automatically contain all the singularities of 
the linear ODEs for $^~ 2) , <" 4) , • • •, ^ 2q) . 

For n = 7, we have checked that the singularities specific to n = 7 (p2 = in 
l|28p . (|29|) ) also appear as singularities of the linear ODE in floating point form (see 
Appendix D for details). For P2 = 1, part of the singularities appear in floating point 

(3) 

form, while for P2 = 2 (i.e. singularities of § H ), no singularities appear in floating 
point form. Similarly, for n = 8, we have checked that the singularities specific to 
n = 8 (p2 = in (|28j) . (|29)l ) also appear as singularities of the linear ODE in floating 
point form (see Appendix D for details). For P2 > 1, no singularities appear in 
floating point form. 

Let us remark that the non observation of some singularities in floating point 

(7) 

form is not really significant. Indeed, we have used 1250 (resp. 1200 terms) for & y H ' 
(resp. <& H ) while the <& H and <& H linear ODEs need more than 20000 terms. 

Figure Q] shows the first family of singularities ([28]) displayed in the complex s 
plane close enough to the unit s-circle. This figure clearly shows a quite rich structure 
for these set of points. This figure looks like a network of nodal points linked together 
by (cardioid-like) curves that can, at first sight, hardly be distinguished from arcs of 
circles. In particular the selected points l + Sw + Aw 2 = as well as the singularities 
for like 1 + 8 w + 20w 2 + 15 w 3 +4w 4 = can be seen to occur quite clearly 

as some of these nodal points. 

f This is a consequence of J23j, J24j, J25J, J26ll yielding k ■ ir = mod. 2tt (see |Appendix E.ij |. 
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Figure 1. First family of singularities l|28|l in the complex s plane (n < 51). 







3- 












1 

o- 

-1- 






-2- 






-3 


-2-10123 



Figure 2. First family of singularities l|28|l in the complex s plane far from the 
unit circle (n < 51). 
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Figure 3. First and second family of singularities J28J, l|29|l in the complex 
plane (n < 16). 



2- 




Figure 4. First and second family of singularities {28j, l|29|l in the complex 
plane far from the unit circle (n < 16). 
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Figure [2] shows the first family of singularities (|28|) far from the unit circle. Figure 
[3] shows all the singularities altogether (first and second family) close to the unit s- 
circle. Finally figure H] shows all the singularities together ( (|28jl . (|29j) ) that are not so 
close to the unit s-circle. 

The accumulation of singularities one can see on figure Q] near s = i and s = — i 
seem to confirm the statement made in Orrick et al [18] that these two points are two 
quite unpleasant points for the susceptibility of the Ising model for which the series 
expansions are not even asymptotically convergent. 

(n) 

Besides reproducing exactly the singularities of the linear ODE for & H , it is 
remarkable to see from the formula l(28|) . |29|) . how to track where each singularity- 
polynomial comes from. This allows one to understand how the singularities of the 
Ising like integrals and (see [7]) and even the Nickelian singularities (J)) 

emerge in these multiple integrals iJH) . This comes simply from the partition l[30|) and 
the equivalent one in 1(28)1 . 



2.2. Singularities: from n-fold integrals to one dimensional integrals 

Consider for instance the singularities 1 — 7w + hw 2 — 4w 3 = occurring in , 
which are given by ((28|) for n = 5, pi = 1 and Vi = 0. As far as conditions on the 
integration angles (see 113311 below), this arises from a situation where two angles are 

(n) 

equal and the three others are equal. Recall that the 3> D integrals are constructed 
with the following restrictions on the angles: 

<t>i = 4>2 = ■■■ = <t>, 4>n = -(n — 1)0. (33) 

One sees that a generalization of this model (f33)) is simply: 

4>i = 4>2 = ■ ■■ = 4>k, 

4>k+i = 4>k+2 = ••• = 4> n , k= 0,1, [n/2]. (34) 

By the condition on the angles, this case is indeed one dimensional, with: 

K = 3 integer. (35) 

The model ((33|) is obviously given by l|34p for k = 1. The Nickelian singularities 
are also given by ((34)) for k = 0, but this time, the underlying model is "zero- 
dimensional". The model constructed along the same lines as in [7] corresponds to 
an integrand: 

n— 1 1 

^ (36) 



3 



The Nickelian singularities arise as poles. 

For k > 2, the singularities given by the model ([34]) . which appear in JH]), are 
thus given neither by ^ nor by Consider one variable of integration such as 

(7]), where the integrand is: 

1 1 

i -x™- 1 ^) -x((n- 1)0) ^ r - ^^^y^^j' ( - 37 - ) 

and denote by such integrals (one then has $^ = ( E > ^^)- 
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Fix n — 5 and k = 2. The constraint ([35)1 on the angles reads: 
3 

<fc = +J 7r ' j integer (38) 

(5) 

with one integration variable. The series of coefficients of <J> 2 is generated along the 
same lines as for (see | Appendix A[ ). The Fuchsian linear differential equation is 
of order six and this order is independent of the value of j in ([38| . The singularities 
of the linear ODE are zeros of the following polynomials: 

w (l- 16w 2 ) (1 + w) (l-3w + w 2 ) (l + 2w-4w 2 ) 

x (l + 4w + 8w 2 ) (l-7w + 5w 2 -Aw 3 ) . (39) 

(5) 

We obtain singularities (from the last two polynomials) appearing for & H and not 

(5) 

occurring for & D . 

The occurrence of the singularities l + 3w + 4w 2 — for (the linear ODE of) $^ 
but not for (the linear ODE of) is explained along similar lines. Note that these 
singularities are common to (the linear ODE of) <P {3) , and . The polynomial 
1 + 3w + Aw appears for (the linear ODE of) from (28]), namely: 

T 2pi (l/2w + 1) = T„_ 2p2 _ 2pi (l/2 W -1). (40) 

One sees that the polynomial 1 + 3w + Aw 2 will appear for all combinations of n, p\ 
and p 2 such that: 

2 Pl =2, n -2p 2 - 2 Pl = 1. (41) 

In other words, the polynomial that arises for given n and pi, will also appear for the 
same value of p\ and for n — 2p 2 . The singularities corresponding to 1 + 3ui + 4w 2 =0 
occur for $^ with n — 5, pi = 1 and p 2 = 1, but lf4lj) is also satisfied for n = 3, 
pi = 1 and P2 = which shows a situation with three angles, with two of them equal. 
This is precisely the integrand in Qj i-e. in . 

Consider now the case n = 6 and k = 2. This amounts to considering the n-fold 
integral with: 

06= -2 0i +J7T, j integer. (42) 

The results are dependent on the integer j. For instance, the series around w = 
reads: 

$( 6) = l + w 6 + 32w 8 ± w 9 + 659 w 10 (43) 

± 1296ui n + 11691 w 12 H 

With the + sign in the series 11431) , the linear differential equation is of order five and 
the singularities are given by the zeros of the polynomials: 

w (l- 16w 2 ) (1 -w) (1 + 2w) (l -9m 2 ) 

x (l-25w 2 ) (l +3w +4w 2 ) . (44) 



The results corresponding to the choice of a minus sign in the series l(43|) are obviously 
obtained b^ w — > — w . We obtain the singularities 1 — 25 w 2 and 1 ± 3u> + Aw 2 = 
occurring for (the linear ODE of) but not for (the linear ODE of) . 

Similarly, for n — 7, (k goes to 3), one obtains for k — 2, the singularities as 
zeros of the following polynomial 1 + 8w + 15w 2 — 21w 3 — 60u> 4 + 16w 5 + 96w 6 + 6Aw 7 , 

§ The last case for n = 6, i.e. k = 3 does not provide singularities other than Nickel's. 
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which has indeed been found numerically in the linear ODE search on a large series 
corresponding to (see | Appendix C[ ). 

We have the remarkable fact that the singularities of the linear ODE for the 

(n) 

multiple integral $ H are given by a finite set of singularities of linear ODEs of 
a set of one-dimensional integrals, namely, N(N + l)/2 one-dimensional integrals, 
with N = [n/2]. For instance, the singularities of the four-dimensional integral <£>^ 
identify with those of, at most, three one-dimensional integrals. This appears, simply, 
from the couple of integers in f28|) which read (2pi, n — 2p2 — 2pi). For fixed n, when 
P2 varies, one sees that we are in fact considering all the lower integer values n — 2pi- 
The same situation holds for |29|) . This identification leads, obviously, to particular 
structures in the singularities for different n. This is what we show in the sequel. 



2.3. Singularity structures of n-fold integrals and particular sets of one- dimensional 
integrals 

The Landau singularities given in | Appendix E are checked against the singularities of 
the linear ODE for <J>^° (n = 3, • • • , 6), and are found to be identical. Assume that 
these formulae do indeed reproduce all the singularities of the linear ODE for , 
for any n. In this case, we can check whether the singularities appearing at n = m 
also occur for n = m + 1, n = m + 2, 

We have found that the singularities at order 2 n will also be singularities at order 
2 n + 2p, where p is a positive integer. Similarly, the singularities at order 2n + l will 
also be present at the following odd orders. 

What is remarkable is the fact that the singularities at odd order also appear at 
even orders. The rule is: all the singularities at odd order n also appear in the higher 
orders (odd and even) except for the first (n — l)/2 even orders. For instance, the 
singularities appearing at n — 3 will occur for all n, except the first even order, i.e. 
4. The singularities appearing at n = 5 will occur for all n, except the first two even 
orders, i.e. 6 and 8. 

The consequence of this embedding of the singularities is the occurrence of some 
singularities at predefined orders. The singularity l + 2w = is present at any order n. 
The singularity 1 — 2w — is present for any even order 2 n. The singularity 1 + w = 
occurs at any order n > 5. The singularity 1 — w = occurs at any order n, except 
for n = 4. All these singularities are Nickelian. The first non Nickelian singularity 
1 + 3w + Aw 2 = appears at all orders n, except for n = 4. 

Moreover, we have given in [7] the Landau singularities for the (linear ODEs 
of the) integrals . These singularities have been found to be identical with the 

(n) 

singularities of the linear ODE for & D obtained exactly up to n = 8 and modulo a 

(n) 

prime up to n = 14. We have seen that all the singularities of the linear ODE of § D 
in the variable s lie in the annulus defined by two concentric circles of radius y/2 and 
The radii of the two concentric circles are the roots, in the variable s, of the 
polynomial 1 + 3w + Aw 2 = 0, that is s 2 + s + 2 = and 1 + s + 2 s 2 = 0. With 
the multiple integrals , one sees that some of the singularities are not confined to 
this annulus anymore. 

Thanks to the Landau conditions, one can now understand this structure from 

(n) (n) 

the reduction of the multiple integrals $ H to a set of one-dimensional integrals $ k 
as far as the location of singularities is concerned. For k = 0, which corresponds to the 
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Nickelian singularities, the "annulus" is the unit circle. For k = 1 corresponding to the 
integrals one has the annulus of radii y/2 and its inverse. For each k, one expects 

the singularities to lie in an annulus with a concentric structure. For these annuli the 
larger radius increases (smaller radius decreases) as k increases. From the reduction of 
the singularities of $^ to these , all the singularities for fixed pi = kin (|28| and 
for fixed p\ — k in ([29)1 will be confined to one annulus. For instance for k = 2, all the 
singularities occurring in the linear ODE for , (i.e. for all n), or, equivalently, all 
the singularities obtained by f28|) for p\ — 1 and by f29|) for p\ — 2 will be confined 
to the annulus of radii 2.79 • • • and its inverse. This value is the root, in the variable 
s, of 1 — 7w + 5w 2 — Aw 3 = occurring for < f > ^' ) . For k — 3, one remarks that the 
annulus will not be obtained from lf28l) which is restricted by 2pi, an even integer. In 
fact this is general. The radii of the annuli are given by ((28j) for k even and by (J29J) 
for k odd. The root in the variable s that will define the annulus occurs at odd order 
n given by 2k + 1. 

The picture now, is as follows. The singularities of the linear ODE for the integrals 
3> H are partionned into "families" indexed by the integer k. The singularities for k = 
are Nickelian and lie on the unit circle, say, r*o = 1. The singularities for k = 1 lie 
in the annulus r% — \/2, l/"\/2 (we discard from now on, the smaller radius). The 
singularities for k — 2 will be confined in the annulus r 2 . The singularities for k = N 
will be in the annulus rjy. These concentric annuli are such that r < r 2 < • • • < r 2 N 
and r\ < r3 < ••• < r 2 N+i, (with r 2 fc < r2k+i)- As fc grows, the radii of two 
neighboring circles behave as r 2 /c+ 2 — f2fc — * and r 2 fe+3 — r 2 fe + i — > 0. This decrease 
is not enough to create an accumulation of circles. We checked with k = 300 circles 
that the decrease goes as k~ a with a < 1 preventing any convergence. For n large 
these radii diverge: rjv — > oo when N — > oo. 

Note that these families, (i.e. the index k) come from the resolution of the Landau 

(n) (ri) 

conditions and from the reduction of the singularities for Q H to the ones of $ k , 
(k = 0, 1, • • • [n/2]). We have no idea as to how these families can be seen directly 
from the mutiple integrals . If the singularities for happen to be identical 
with those occurring in the linear ODE for x {n) , it may become important to see 
whether this picture persists and whether this picture is showing another partition of 
the susceptibility \ instead of the known sum on \^ ■ 

Figure [U [6] and [7] show how the first family of singularities (|28| in the s complex 
plane is decomposed according to the integer k in ([34]) . Figure 5 shows singularities 
(|28p for a given odd value of k, namely k = 5 for any odd values of n up to 91. 
Figure [6] shows singularities lf2"8| for a given even value of k, namely k = 2 for any 
odd values of n up to 71. Figure 7 shows singularities f28|) for a given even value 
of k, namely k — 6 for any even values of n up to 80. The figures corresponding 
to the filtration of the singularities of the first family (J28J) in terms of the integer k 
(previously displayed altogether with figures Q] and [2]) deserve some comments. First, 
one sees that the various "crescent" corresponding to different values of k are very 
similar. Secondly one sees from figure 5 that the odd n, odd k "crescent" break the 
s <-> — s symmetry (for even n, even k, the equations for the set of singularities are 
functions of s 2 , see figure[7| in a quite dramatic way: the singularities in the "crescent" 
of figure [5] all lie only in the left half s-complex plane. Similarly the singularities in 
the "crescent" of figure [6] all lie in the right half s-complex plane. 

Along this s «-> —s symmetry line it is worth recalling that the low-temperature 
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Figure 5. Crescent in the complex s plane given by l|28ll : k = 5, n < 91, n 
odd. 
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Figure 6. Crescent in the complex s plane given by l)28ll : k = 2, n < 71, n 
odd. 
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Figure 7. Crescent in the complex s plane given by l|28ll : k = 6, n < 80, n 
even. 

susceptibility of the Ising model has this s <-> — s symmetry (the low-temperature 
susceptibility is a function of s 2 or w 2 ) but the high-temperature susceptibility breaks 
that s «-* — s symmetry, and this is also the case for the rt-fold integral with 
n odd. Our rt-fold integrals are introduced to provide an educated guess as to 
the location of the singularities of the X ■ As far as location of singularities of the 
are concerned, it is not totally clear for n odd if the s «-> — s (resp. w <-> — w) 
symmetry will not be partially restored on the global set of singularities with the 
occurence for a singularity P n (w) = for a given value of n, of the opposite value 
for, perhaps, a different value of n : P m (—w) = 0. 

Remark: Quite often, in this paper, we use (by abuse of language) the words 
"singularities of an n-fold integral" to describe a larger set of singularities, namely 
the singularities of the linear ODEs that the rt-fold integral satisfies. A rigorous 
study would require, for any "singularity", to perform the (differential Galois group 
and connection matrix) analysis we have performed in [6]. It amounts to getting 
extremely large series, deduced from the obtained linear ODE, that coincide with 
the series expansion of the n-fold integral we are interested in, and find out if these 
series actually exhibit these singularities. With this tedious, but straightforward, 
procedure we can extract the singularities of a specific n-fold integral among the 
larger set of singularities of the corresponding linear differential equation. In view of 
the large number of singularities we display in this paper, we have not performed such 
a systematic analysis, that would have been quite huge. Furthermore it is important to 
note that this "connection matrix" approach [6] requires to have the linear ODE of the 
n-fold integral. A knowledge of the linear ODE modulo a prime is not sufficient. We 
could have performed this analysis for $^ and , but, in that case, we already 
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have a deeper result [6] namely the connection matrix analysis for x an d 
providing an understanding of the singularities of these n-fold integrals themselves 
(in w and also in s). 

Right now, the only singularities found for the x > other than Nickelian, are 
the quadratic roots of 1 + 3w + Aw 2 = 0, (i.e. the first annulus) which appear at all 
orders (except n = 4 for Let us show, in the sequel, how this polynomial can 

be "special". 



3. Towards a mathematical interpretation of the singularities 

In a set of papers [19} [20] , we have underlined the central role played by the elliptic 
parametrization of the Ising model, in particular the role played by the second order 
linear differential operators corresponding to the complete elliptic integral E (or K), 
and the occurrence of an infinite number of modular curves [12j, canonically associated 
with elliptic curves. The deep link between the theory of elliptic curves and the theory 
of modular forms is now well established [2 lj . 

Consequently, it may be interesting to seek "special values" of the modulus 
k, (singularities of the X ) that might have a "physical meaning", as well as a 
"mathematical interpretation". 

For that purpose, recall that the modular group requires one to introduce the 
elliptic nome, defined in terms of the periods of the elliptic functions, 

/ K{l-k 2 )\ f . , 

q= exp^-7r K , k 2\ ) = ex P( i7rr ) (45) 

and the half period ratic0 r. We write the complete elliptic integral K as 

K(k) = 2 F 1 (l/2, 1/2; l;fc). (46) 

Relations between K(k) evaluated at two different moduli can be found in, e.g. [22J. 



3.1. Some isogenics of elliptic curves seen as generators of the renormalization group 

The arguments in K in these identities are related by the so-called, respectively, 
descending Landen and ascending Landen (or Gauss) transformations: 



= 1 - VT ^ (47) 

* - h - ^ (4«) 

These transformations (or correspondences [23] [24]). decrease or increase the modulus 
respectively. Iterating lj47j) or lf48]) . one converges to k = or k = 1 respectively. 
The half period ratio transforms through ([47]) , ([48]) , as 



t^2t, t -> -t (49) 

respectively. The real fixed points of the transformations (|47| and ([48]) are k = 
(the trivial infinite or zero temperature points) and k = 1 (the ferromagnetic and 

If In the theory of modular forms q 2 is also sometimes used instead of q. In number theory literature 
the half-period ratio is taken as —it. 
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antiferromagnetic critical point of the square Ising model) . In terms of the half period 
ratio, this reads r = oo and r = respectively, which correspond to a degeneration 
of the elliptic parametrization into a rational parametrization. In view of these fixed 
points, it is natural to identify the transformations l(47j) or (|48|) . and more generally 
any transformation^ r — > n ■ r or r — » r/n (n integer), as exact generators of the 
renormalization group of the two-dimensional Ising modetfl. 

One does not need to restrict the analysis to the real fixed points of the 
transformations. If one considers the Landen transformation (|38)l as an algebraic 
transformation of the complex variable k and if one solves k\ — k 2 = 0, one obtains: 

fc-(l-fc)-(fc 2 + 3fc + 4) = 0. (50) 

The quadratic roots 

k 2 + 3k + A = 0, (51) 

are (up to a sign) fixed points of (|48|) . We thus see the occurrence of additional 
non-trivial complex selected values of the modulus k, beyond the well-known values 
k = 1, 0, oo (corresponding to degeneration of the elliptic curve into a rational curve). 
Physically, these well-known values k = 1, 0, oo correspond to the critical Ising model 
(k — 1) and to (high- low temperature) trivializations of the model (A; = 0, oo). 

3.2. Complex multiplication for elliptic curves as (complex) fixed points of the 
renormalization group 

We come now to our point. The first "unexpected" singularities 1 + 3 w + 4 w 2 = 
found [4j [5] for the Fuchsian linear differential equation of > an d also in other 
n-fold integrals of the Ising class [7] , reads in the variable k = s 2 as 

(fc 2 + 3k +4) (4k 2 + 3k +1) = 0. (52) 

The first polynomial corresponds to fixed points of the Landen transformation (see 
l|50p). In other words we see that the selected quadratic values I + 3w + 4w 2 = 0, 
occurring in the (high-temperature) susceptibility of the Ising model as singularities of 
the three-particle term x' 3 ^ , can be seen as fixed points of the renormalization group 
when extended to complex values of the modulus k. 

For elliptic curves in fields of characteristic zero, the only well-known selected 
set of values for k corresponds to the values for which the elliptic curve has complex 
multiplication |26J . Complex multiplication for elliptic curves corresponds to algebraic 
integer values (integers in the case of the Heegner numbers, see |Appendix F[ ) of the 
modular j-function, which corresponds to Klein's absolute invariant multiplied by 
(12) 3 = 1728: 

(l-k 2 + k 4 ) 3 , x 

j(k) = 256 --^ (53) 

1 ^ fc 4 -(l-fc 2 ) 2 V ; 

A straightforward calculation of the elliptic nome lf45j) gives, for the polynomials 
l|52p . respectively, an exact value for r, the half period ratio, as very simple quadratic 
numbers: 

±3+iV7 ±l+iV7 
n = ^ , t 2 = (54) 

§ See relation (1.3) in [25] . 

\ A similar identification of these isogenies r — > n ■ r with exact generators of the renormalization 
group can be introduced for any lattice model with an elliptic parametrization (Baxter model, ...). 
| Note that the two polynomials in l|52|l are related by the Kramers- Wannier duality k — > l/k. 
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These quadratic numbers actually correspond to complex multiplication of the elliptic 
curve and for both one has j = (— 15) 3 . These two quadratic numbers are such that 
2 n =F 1 = T 2- Let us focus on t 2 for which we can write: 

2 , s 

r = 1 - -. (55) 

T 

Taking into account the two modular group involutions r — > 1 — r and 
t — > 1/r, we find that 1 — 2/r is, up to the modular group, equivalent to r/2. 
The quadratic relation t 2 — r + 2 = thus amounts to looking at the fixed points of 
the Landen transformation t ^ 2t up to the modular group. This is, in fact a quite 
general statement. The complex multiplication values can all be seen as fixed points, 
up to the modular group, of the generalizations of Landen transformation, namely 
t — ► n t for n integer, t 2 — t + n = or t = 1 — ~ rt • r, where ~ denotes the 
equivalence up to the modular group. 

|Appendix G presents an alternative view by considering the solutions as fixed 
points under Landen transformations of the modular j— function. 

In view of the remarkable mathematical (and physical) interpretation of the 
quadratic values l + 3w + 4w 2 = Oin terms of complex multiplication for elliptic 
curves, or fixed points of the renormalization group, it is natural to see if such 
a "complex multiplication of elliptic curves" interpretation also exists for other 
singularities of x > and as a first step, for the singularities of the linear differential 
equations of our n-fold integrals iJH), that we expect to be identical, or at least have 
some overlap, with the singularities of the X • 

Noting that the modular j— function is a function of s 2 or w 2 (see (|F,2|) in 
|Appcndix F [ the occurence of l+3w + 4w 2 = Oasa selected quadratic polynomial 
condition means, at the same time, the occurence of the other quadratic polynomial 
condition l-3w +4w 2 = (see |Appendix F| and |Appendix G.2[ ) . 

Besides l-3w+4w 2 = 0, we have found two other polynomial conditions which 
correspond to remarkable integer values of the modular j-function. The singularities 
1 — 8w 2 = correspond to j = (12) 3 and r = ±1 + i (see | Appendix Ft . They 
correspond to "Nickelian singularities" for x*" 8 ^ (and thus ) and to "non-Nickelian 
singularities" for < &^°' ) . Another polynomial condition is 1 — 32 w 2 — , which 
gives "non-Nickelian singularities" that begin to appear at n — 10 for < £ > ^°' ) - These 
singularities correspond to the integer value of the modular j-function, j = (66) 3 and 
to r = 2 % or t == -4/5 + 2i/5. 



3.3. Beyond elliptic curves 

Among the singularities of the linear ODE for $^ given in lfT2|) . lfT3|) or obtained from 
the formula given in |Appendix E| up to n — 15, we have found no other singularity 
identified with selected algebraic values of the modular j-function corresponding 
to complex multiplication for elliptic curves. Could it be that the (non-Nickelian) 
singularities lfT2]). <fl3)) . which do not match with complex multiplication for elliptic 
curves, are actually remarkable selected situations for mathematical structures more 
complex than elliptic curves ? With these new singularities, are we possibly exploring 
some remarkable "selected situations" of some moduli space of curves corresponding 
to pointed (marked) curves [27], instead of simple elliptic curves [28]? In practice 
this just corresponds to considering a product of n times a rational, or elliptic, curve 
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minus some sets of remarkable codimension-one algebraic varieties [TT], x% Xj = 1, 
XiXjXk — 1, hyperplanes Xi — xj , ■ ■ ■ 

We try to fully understand the singularities of the n-fold integrals corresponding 
to the x i that is to say particular n-fold integrals linked to the theory of elliptic 
curves. These n-fold integrals are more involved than the (simpler) n-fold integrals 
introduced by Beukers, Vasilyev (29] [30] and Sorokin [31] [32], or the Goncharov-Manin 
integrals |33j which occur in some moduli space of curves [34, 35j simply corresponding 
to a product of rational curves (CP\ xCPi ■ ■ ■ xCP\). An example of such integrals, 
linked^ to £(3), is displayed^ in | Appendix H] 

Let us close this section by noting that Heegner numbers and, more generally, 
complex multiplication have already occurred in other contexts, even if the statement 
was not explicit. In the framework of the construction of Liouville field theory, Gervais 
and Neveu have suggested [UJ new classes of critical statistical models, where, besides 
the well-known A^-th root of unity situation, they found the following selected values 
of the multiplicative crossing t [42] : 

t = ei .(iW3)/2 = (56) 
t = e l77{1+l) = -e"\ (57) 

If one wants to see this multiplicative crossing as a modular nome, the two previous 
situations actually correspond to selected values of the modular j-function namely 
j((l + i%/3)/2) = (0) 3 for (56]), and j(l + i) = (12) 3 for (5Z|), which actually 
correspond to Heegner numbers and, more generally, complex multiplication |26j . It is 
however important not to feed the confusion already too prevalant in the literature, 
between a "temperature-like" nome such as (|45ll and a multiplicative crossing modular 
nome. In the Baxter model [43] [44], the first is denoted by q and the second one 
by x. In fact one probably has, not one, but two modular groups taking place, one 
acting on the "temperature-like" nome q and the other acting on the multiplicative 
crossing x. We will not go further along this quite speculative line which amounts to 
introducing elliptic quantum groups [45] and elliptic gamma function^ (generalization 
of theta function£tj]) • 



4. Conclusion 

The ultimate goal of our "Ising class" integrals is to get some insight into the x^ n ' 
and, hopefully, into the susceptibility of the Ising model. For that purpose we 
have introduced n-fold integrals J8]) such that we expect the singularities of the 
corresponding linear ODE to overlap, as much as possible, with the singularities of 
the linear ODE for the x ■ We have obtained the linear differential equations for 

% Note that C ( or the polyzeta) function evaluated at integer values (C(3), C(5), ■■■) do occur in our 
more involved n-fold integrals, in particular in the representation of the connection matrices [5] of 
the differential Galois group of the Fuchsian linear ODEs of x ■ 

tt These n-fold integrals 06, .37, 38, 39, go] look almost the same as the ones we have introduced 
and analyzed in the study of the diagonal susceptibility of the Ising model for which n-th root 
of unity singularities occur. 

f Which can be seen [46] as "automorphic forms of degree 1" when the Jacobi modular forms are 
"automorphic forms of degree 0" and are associated (up to simple semi-direct products) with SL(3, Z) 
instead of SL(2, Z) 

ft The partition function of the Baxter model can be seen as a ratio and product of elliptic gamma 
functions and theta functions. It is thus naturally expressed as a double infinite product. Similar 
double, and even triple, products appear in correlation functions of the eight vertex model [47 48J. 
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these n-fold integrals <& H , up to n = 4 and up to n — 6 modulo a prime. Prom these 
exact results together with an exhaustive Landau singularity analysis, we provided a 
quite complete description of the singularities of these linear ODEs. 

From the Landau conditions, the singularity structures are explained. The 
singularities corresponding to $ H are found to also occur at a higher predefined order 
p > n. With these multiple integrals and the associated Landau conditions, we have 

(n) 

been able to understand why the simple integrals $ D have succeeded reproducing the 
Nickelian singularities and the new quadratic l+3w+4w 2 — 0. These simple integrals 

(n) 

appear to be "a first approximation" to & H . Other one-dimensional integrals pop up 
to account for the additional singularities not occurring for . 

We have then a remarkable finding that, the singularities for the multiple integrals 
can be associated with the singularities for a finite number of one dimensional integrals. 
If the singularities, associated with these n-fold integrals ([8]), happen to be identical 
with (or to overlap) the singularities associated with the x^ > it becomes important to 
understand this mechanism for the x themselves. If this mechanism of singularity 
embedding occurs for x > it might be explained by a Russian doll structure for the 
same linear differential operators. We know that the linear differential operator for 
X- 1 ' (respectively X ) is "contained" in (rightdivides) the linear differential operator 
for x^ (respectively x^ 4 -*), and furthermore we even have direct sum decomposition 

(n) 

properties. For the 0> H , it is not these mechanisms which are at work. 

Our primary goal in this study is to identify as many singularities as possible 
for the x^ • The singularities of the ODEs associated with the $^ quantities 
correspond, in the Landau equations framework, to leading pinch singularities 
(relatively to the singularities manifolds D(<f>, ip) = 0). For the other quantities 
previously studied [7] which belong to the Ising class integrals, the same feature holds. 

At this step, the natural questions arising are: whether the scheme, from the 
Landau singularities point of view, which holds for still holds for x^ an d 

whether the singularities of can be considered as singularities of the x^ ? 

From the Landau singularities viewpoint, the Fermionic determinant G(n) 2 is 
going to introduce new manifolds of singularities. When the Lagrange multipliers 
relative to the singularities manifolds introduced by the Fermionic determinant are all 

(n) 

set equal to zero, one deals with the Landau equations of the & H quantities. Thus 
the singularities obtained for the quantities are also solutions of the Landau 

equations of the X ■ However this feature does not mean that the singularities of 
the <I>^ quantities will necessarly appear as singularities of the x^ ODEs. Indeed 
some selection rules may take place and may reject some of them. For instance, one 
expects singularities linked to the Y[ Hi to occur for the Landau singularities of the 

(n) 

$ H . One finds that some selection rules exclude them. Our "educated guess" is 

(n) 

that all the Landau singularities of the $ H will be in the Landau singularities of 
the X j however we do not exclude the possibility that the x will have more 

(n) 

Landau singularities than the & H . Another "educated guess" is that the Landau 
singularities of the x will exhibit a similar embedding that the one we found for 
the ^\ This naturally raises the question already considered in [8], of a "strong" 
Russian doll structure for the linear differential operators of the x^ > namely that the 
linear differential operator of (resp. x' 4 ^) could right-divide the linear differential 
operator of x^ (resp. x^)j an d so on. 
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This knowledge of the singularities will help in the search for the corresponding 
linear ODE. For instance, we have 24 head polynomial "candidates" for x' 5 ' an d 19 
"candidates" for x that can, from the outset, be put in front of the higher order 
derivative of the unknown linear ODE. From the knowledge we have gained from all 
these n-fold integrals of the "Ising class", one can guess the order of magnitude of the 

(5) 

multiplicity of some singularities. Furthermore, as shown for the linear ODE for 
and (and also from previous ODEs), we know that the "cost" (in terms of the 

number of series coefficients) will be much less for a non minimal order linear ODE 
than for the minimal order one. 

(n) 

Concerning the non Nickelian singularities that the multiple integrals & H have 
given, we focussed on 1 + 3w + 4w 2 = which actually occurs for the linear ODE 
of x' 3 -*, or for x^ seen as a function of s. As far as a mathematical interpretation is 
concerned, we have shown that this quadratic polynomial condition corresponds to a 
selected situation for elliptic curves namely the occurrence of complex multiplication. 
The other non-Nickelian (candidate) singularities, (fl"2|) . lfl3|) do not correspond to 
complex multiplication of elliptic curves. 

Assuming that the non Nickelian singularities obtained in the linear ODE for the 
integrals J8]), will be, at least, included in those for the X i various lines of thought 
are possible. 

One may imagine that the decomposition of the susceptibility of the Ising model 
in terms of an infinite sum of x ls quite an artificial one with no deep mathematical 
meaning, i.e. x are quite arbitrary n-fold integrals. In this case, no interpretation 
within the theory of elliptic curves has to be looked for and the occurrence for 
1 + 3w + 4w 2 = Oof complex multiplication for elliptic curves would be just a 
coincidence. 

Another option amounts to saying that one needs to introduce (motivic) 
mathematical structures [36, 37, 38, 39, 40j beyond the theory of elliptic curves (moduli 
spaces, marked curves, ...), and beyond the elliptic curves of the Ising (or Baxter) 
model, to get a mathematical interpretation of these singularities. We tend to favour 
the latter option. 
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(n) ■ . ^.(n) 

Appendix A. Series expansions o{<P H and of single integrals <Pj, 

We give in this Appendix, the series expansion that has been used for <I>^ . Expanding 
the integrand of in the variables Xj , one obtains 

1 n-1 r 2Tr 11 oo n 
j=l J p= o i=l 

We make use of the yj x? Fourier expansion [H [8] 

oo 

y jX P = £ J*l . a (p,|fc|)-4\ Z^expfi^-) (A.2) 

/c— — oo 

where a(k,p) is a non-terminating hypergeometric function that reads (with to = 
a(k,p) = ( T f) x (A.3) 



We define (p) by 



• 1 + m 1 + to 2 + m 2 + to 2 
4-F3 — — > — 7 — , — 7 — , — 7 — ; 1 + k, 1 + p, 1 +to; 16 w 



where the angular constraint is introduced through the delta function that has the 
Fourier expansion: 




= J2 {ZiZ2---Z n ) k (A.5) 

/ k— — oo 

The integrals (|A.1|) become 

1 00 00 n 

-• Y, £(2-^o)-<II^4> ( A -6) 



00 00 

H Til 

k—~oo p—0 3=1 



where the integration is over independent angles. 

Using the Fourier expansion (| A.2[) . one obtains the integration rule 

(yjtfZt) = w^-a(p,\k\) (A.7) 

and finally: 



00 00 

H n. 

k=0 p=0 



00 00 

rEE( 2 -^»)'(2-^,o)-^ +p) -a n (k,p). (Ai 



The derivation of the series expansions for the one dimensional integrals ((37] 
proceeds along similar lines. The integrand of the integrals IpT)) is expanded in x 



with 



n — k 27rj 
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Here, we use the Fourier expansion 

oo 

x m = w m ■^2(2-5 p ,o)-w p ■b(p,m)-cos{p(t>) (A.ll) 

where b(k,p) is a non-terminating hypergeometric function that reads (with m = 
k+p): 

b(k,p) = ( m 7 l \ x (A.12) 



4-T3 



l + ml + m2 + mm 2 

■ 1 + re, 1 -\- p, 1 + rn; 16 w 



2 7 2 ' 2 ' 
The integration of the one-dimensional integrals (|A.9|1 gives 



oo oo oo 



with 



and 



where 



v ' u ' p=0 pi=0 P2=0 

x wP»+pi+» .b(jn,p{n-k)) -b(p 2 ,pk) I(pi,p 2 ) (A.13) 



I(pi,P2) = 2 ' ( 1 + V.o)-cos(c), for p 2 -(n-k)= k-pi, 



1 6 2 

I{Pi,P2) = — ^ ' sin(&7r) • cos(&7r — c), for p 2 ■ (n — k) ^ k ■ px, 



pi 



b=—^-p 2 , c= ~j L -P2- (A.14) 

Appendix B. Linear differential equations of some <f>^ 

(3) 

Appendix B.l. Linear ODE for <& H 

(3) 

The minimal order linear differential equation satisfied by & H reads 



where 



with 



dw 

n=0 



a 5 (w) = (l-w)(l-4w) 4 (l + 4w) 2 (l + 2 l «) 

x (l + 3w + 4w 2 ) -w 3 -P 5 H, (B.2) 
a 4 (w) = (1 - 4io) 3 (1 + 4w) ■ w 2 ■ P 4 (w), 

03(10) = -2 (1 - 4w) 2 • iuP 3 (iy), a 2 (iu) = 2 (1 - Aw) ■ P 2 (w), 
ai (w) = -8P X H, a H = -96P H, 

P 5 (w) = -5 + 21 w + 428 w 2 + 5364 w 3 - 82416 w 4 - 299504 
+ 714944 w 6 + 3127872 w 7 - 8220672 w s - 25858048 w 9 
- 7077888 w 10 + 31424512 w 11 - 42467328 w 12 
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- 31457280 w 13 - 4194304 w 14 + 4194304 w 15 , 

P 4 (w) = -40 + 7 + 5232 w 2 + 37159 w 3 - 447778 w 4 - 4947500 w 5 
+ 19493448 w 6 + 258464112 w 7 + 499205984 w 8 - 1612751808 w 9 

- 4667817856 w 10 + 13827459072 w 11 + 67078416384 w 12 

+ 62392041472 w 13 - 81535369216 w 14 - 116835483648 w 15 
+ 124662054912 w 16 + 146016305152 w 17 - 197258117120 w 18 

- 131667591168 w 19 - 11676942336 w 20 + 150 3 2 3 8 5 5 36 w 21 , 
P 3 (w) = 35 - 25 w- 8683 w 2 - 10149 w 3 + 619246 w 4 + 5273820 w 5 

- 52472072 w 6 - 588147792 w 7 + 491073248 w s + 18721819584 w 9 
+ 47622771584 w 10 - 97459630592 w 11 - 441418588160 w 12 

+ 651003559936 w 13 + 4694018588672 w 1A + 4729946636288 w 15 

- 7193770328064 w 16 - 11814519701504 w 17 + 7399599505408 w 18 

+ 10981996494848 w 19 - 16439524196352 w 20 - 10434623045632 w 21 

- 916975517696 w 22 + 1125 28 1 43 1 552 w 23 , 

P 2 (w) = -10 + 101 w + 11088 w 2 - 42855 w 3 - 1117278 w 4 - 1918516 w 5 
+ 72221464 w 6 + 460080656 w 7 - 4999186016 w s 

- 33474428224 w 9 + 67440200320 w 10 + 808560558592 w 11 

+ 535166693376 w 12 - 6771457933312 w 13 - 7468556451840 w 14 
+ 46143514476544 w 15 + 91488863125504 w 16 

- 75107733078016 w 17 - 239438663778304 w 18 + 31904728350720 w 19 

+ 234058806198272 w 20 - 237446193217536 w 21 - 164181567340544 w 22 

- 18975165513728 w 23 + 16973710753792 w 24 , 

Pi(io) = -5- 1142 w + 8106 w 2 + 210846 w 3 - 1070376 w 4 - 7771160 w 5 

- 22029952 w 6 + 833894752 w 7 + 3334510976 w 8 - 39736449920 w 9 

- 156101859328 w 10 + 663306718208 w 11 + 2995615555584 w 12 

- 5033154314240 w 13 - 26250785980416 w 14 + 28618066755584 w 15 

+ 158047775227904 w 16 - 42836217036800 w 17 - 410317620248576 w 18 

- 95925074657280 w 19 + 462245318361088 w 20 - 328990199906304 w 21 

- 249443110617088 w 22 - 35270271434752 w 23 + 24464133718016 w 24 , 
P (w) = -5 + 58 w + 3234 w 2 - 18994 w 3 - 229330 w 4 + 1516 w 5 

+ 7017504 w e + 74689472 w 7 - 647069792 w 8 - 4260373952 w 9 
+ 15887163648 w 10 + 96789618688 w 11 - 136120508416 w 12 

- 917765144576 w 13 + 877996605440 w 14 + 5646695006208 w 15 

- 2888887697408 w 16 - 16785155817472 w 17 - 5241017729024 w 18 
+ 17952426426368 w 19 - 13058311192576 w 20 - 9329742708736 w 21 

- 1275605286912 w 22 + 824633720832 w 23 . (B.3) 
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Appendix B.2. Linear ODE for 

The minimal order linear differential equation satisfied by reads (with x — 16w 2 ) 
6 rP 

E a ^ ■ dx^ F W = °> ( B - 4 ) 

n=0 



where 



with: 



a 6 {x) = 64 (x - 4) (1 - a;) 4 a; 4 • P 6 (x), a 5 (x) = -128 (1 - x) 3 x 3 ■ P 5 (x), 
a 4 {x) = 16 (1 - x) 2 x 2 ■ P 4 (x), a 3 (x) = -64 (l-x)x- P 3 (x), 
a 2 {x) = -A- P 2 {x), ai(x) = -8-P^x), a Q (x) = -3 (1 - x) ■ P {x), 

P 6 (x) = 128 + 2233 x - 2847 x 2 + 3143 x 3 - 3601 x 4 + 144 a; 5 - 64 a; 6 , 
P 5 {x) = 3712 + 51523 x - 216377 a; 2 + 289918a; 3 - 312896a; 4 

+ 262111 a; 5 - 63167a; 6 + 5512 x 7 - 896 a; 8 , 
Pi(x) = -121856- 1102304a; + 11038289a; 2 - 26106487a; 3 

+ 31515802 x 4 - 31027694 x 5 + 21291429 x 6 - 5166011 x 7 

+ 410160 a; 8 - 67776a; 9 , 
P 3 (x) = 38144 + 10604 x - 4644281 a; 2 + 20909702 a; 3 - 37890772 x 4 

+ 42011874 a; 5 - 37552559 x e + 22474036 x 7 - 5465572 a; 8 

+ 392536 a; 9 - 65984 a; 10 , 
P 2 (x) = 163840 - 4162688 a;- 18120152 a; 2 + 277110610 a; 3 

- 880048289 x 4 + 1357147519 a; 5 - 1395938590 x e + 1 141353668 x 7 

- 621323833 a; 8 + 150842795 a; 9 - 9676720 a; 10 + 1656512 a; 11 , 
Pi (a;) = -366592 + 3113752 x + 17465700 a; 2 - 120658444 a; 3 

+ 240321805 a; 4 - 259277988 a; 5 + 219951814 a; 6 - 142314304 a; 7 
+ 42534921 a; 8 - 2056040 a; 9 + 435200 a; 10 , 
P (x) = 561152- 1496400a;- 13171575 a; 2 + 30840556 a; 3 - 24381 198 a; 4 
+ 20352948 a; 5 - 13268091 a; 6 + 309360 x 7 - 120000 a; 8 . 



Appendix B.3. Linear ODE modulo a prime for $ H 

The 1 
form 



(5) 

The linear differential equation of minimal order seventeen satisfied by $> H is of the 



n=0 



with 



a 17 (w) = (l-4w) 12 (l + 4w) 9 (l-w) 2 (u; + l)(l + 2w) 
x (l + 3w + Aw 2 ) 2 (l - 3w + w 2 ) (l + 2w-4w 2 ) 
x (1 +4w + 8w 2 ) (1 - 7w + 5w 2 - 4w 3 ) 



x 



(l-w-3w 2 +4:w 3 ) (l + 8w + 20w 2 + 15w 3 + 4w 4 ) • w 12 • P 17 (w), 
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a 16 (w) = w 11 (1 - 4w) n (1 + 4w) 8 (1 - w) (l + 3w + Aw 2 ) ■ P 16 (w), 
a 1B (w) = w 10 (l-4to) 10 (l + 4j t i) 7 -P 15 H 
ou (to) = w 9 (1 - 4w) 9 (1 + 4wf ■ P m (w), 

where the 430 roots of Pn(w) are apparent singularities. The degrees of 
these polynomials P n (w) are such that the degrees of a, (to) are decreasing as: 
deg(di+i(w)) — deg(ai{w)) + 1. In fact, with 2208 terms we have found the ODE 
of $^ 5) at order q = 28 using the following ansatz for the linear ODE search (Dw 
denotes d/dw) 

Q 

»(*) • p(») • Dw 1 (B.6) 

*=o 

with: 

a (i) = ^(-W) • (1 - 16 W 2 ) Q (" 1+i ) • s ^ (1+J - 9) (B.7) 
where a(n) = Mm(0, n) and 

s = (1 + to) • (1 - w) ■ (1 + 2 w) ■ (1 - 3 w + w 2 ) (1 + 2 u> - 4 w 2 ) 

x (1 + 3w + Aw 2 ) ■ (1 + Aw + 8w 2 ) ■ (1 - 7w + 5w 2 - Aw 3 ) 

x (1 - w - 3 w 2 + 4 w 3 ) • (1 + 8 w + 20 w 2 + 15 w 3 + 4 w 4 ) 

the p(i) being the unknown polynomials. 

The minimal order ODE is deduced from the set of linear independant ODEs 
found at order 28. 

Appendix B.4- Linear ODE modulo a prime for *g> 

The linear differential equation of minimal order twenty-seven satisfied by $^ rea ds 
(with x = w 2 ) 

E «»(*) • - o, (B.s) 

n=0 

with 

027(2;) = (l-16x) 16 (l-4x) 3 (l-x)(l-25x)(l-9x)x 21 

x (l -x+ 16a; 2 ) (l - 10x + 29x 2 ) • P 2 7(x), 
M*) = (l-16x) 15 (l-4x) 2 x 20 • P 26 (x), 
a 25 (x) = (1 - 16 x) 14 (1 - 4a;) x 19 • P 25 (z), 
o 24 (x) = (1- 16 x) 13 x 18 • P u (x), 

(B.9) 

where the 307 roots of P 2 7(x) are apparent singularities. The degrees of the P n {w) 
polynomials are such that the degrees of a* (to) are decreasing as: deg(oj+i(to)) = 
deg(a,i(w)) + 1 

In fact, with 1838 terms we have found the linear ODE of at order q = 42 
using the following ansatz for the linear ODE search (Dx denotes d/dx) 

Q 

s{i) ■ p(i) ■ Dx l (B.10) 
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with: 



s(i) = x a <-- 1+ ^ ■ (1 - 16x) a <-- 1+ ^ ■ s^ 1+l - q) 



(B.11) 



where a(n) 



Min(0, n) and 

s = (1 - 25 z) • (1 - 9.x) • (1 -Ax) ■ (1 - x) 



x (l-a; + 16a; 2 )-(l-10a; + 29a; 2 ) 



(B.12) 



the p(i) being the unknown polynomials. 

The minimal order ODE is deduced from the set of linear independant ODEs 
found at order 42. 

(7) (81 

Appendix C. Singularities in the linear ODE for & H and & H 

(7) 

For & H , we generated large series, (1250 coefficients and 20000 coefficients modulo 
primes), unfortunately, insufficient to obtain the corresponding linear ODE. However, 
by steadily increasing the order q of the ODE (and consequently decreasing the degrees 
n of the polynomials in front of the derivatives), one may recognize, in floating point 
form, the singularities of the ODE as the roots of the polynomial in front of the higher 
derivative. A root is considered as singularity of the still unknown linear ODE, when 
as q increase (and consequently decreasing n), it persists with more stabilized digits. 

(7) 

Using 1250 terms in the series for & H , the following singularities are recognized 
(1 -4 to) (l - 5w + 6w 2 - w 3 ) (l + 2w-8w 2 -8ic 3 ) (1 + Aw) ■ w 



(l + 2w-w 2 -w 3 ) (l-3w + w 2 ) (l + 2w-4w 2 ) (1 + w) 

(16 w s - 32 w 7 - 17 w 6 + 62 w 5 - 5 w 4 - 35 w 3 + 10 w 2 + 3 w - l) 

(64 w 7 + 96 w 6 + 16 w 5 - 60 w 4 - 21 w 3 + 15 w 2 + 8 w + l) 

(128 w 5 - 32 w 4 + 48 w 3 + 16 w 2 + 4 w - l) 

(4 w 5 + 51 w 3 - 21 w 4 - 1 + 10 w - 35 w 2 ) 

(Aw 3 + 7w - 5 w 2 - 1) (Aw 4 + 1 + 7w + 26w 2 + 7w 3 ) 

(A w 4 + 1 + 8 w + 20 w 2 + 15 w 3 ) 

(l + 12w + 54w 2 + 112w 3 + 105w 4 + 35w 5 + 4w 6 ) = 0. 



Note that we have not seen with the precision of these calculations the occurence 

(3) 

of the singularities of the & H . 

(8) 

With similar calculations using 1200 terms for ^ , the following singularities 
are recognized: 



(l-2w)(l + 2w) (l - 2u> 2 ) (1 - Aw) (l - Aw + 2w 2 ) (1 + Aw) 
(l + Aw + 2w 2 ) (8w 2 - l) (3w- 1) (1 - w) (1 + w) (3w+ l)w 
(1138 w w - 1685 w s + 960 w e - 242 w 4 + 26 w 2 - l) 
(32 w 4 - 10 w 2 + 1) (1312 w e - 56 w 4 + 30 w 2 - l) 
(10 w 2 - 6w + 1) (4w 3 - 8w 2 + 6w- l) 
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(5w-l)(l+2w 2 ) (5w + l) 

(l0u; 2 + 6w + l) (4w 3 + 8w 2 + 6w+ 1) = 0. 

We will see in the next section [Appendix E.3| that we missed the polynomials: 
(l-3w + 4w 2 ) (l + 3w + 4w 2 ) (l-10w 2 + 29w 4 ) 

Note that the stabilized digits in these singularities can be as low as two digits. 

(n) 

Appendix D. Landau conditions and pinch singularities for <f> D and 
integrals of LJ Vj • 

Similarly to the integral representation (fl6| of one has: 

[ 2 * def, dj> y/(l - 2w COS 0)2 - 4^ 

Jo 2* J o Jfa D(<f>^) (DA) 

- 2w cos((n - 1)0)) 2 - 4w 2 

X 



and 



D((n-1)4>, (n-l)V) 



For the singularities of the associated ODEs are given as solutions of: 

D{4>, 1>) = 0, 

D{{n-1)4>, (n-l)^) = 0, 

ai • sin(0) + 0:2 • sin((n — 1) 0) = 0, with a\, a-i ^ 0, 

ai ■ sin(V') + a 2 • sin((n - 1) ^) = (D.3) 

which are nothing less than the Landau conditions restricted to pinch singularities 
of the singularity manifolds D{4>i, ipi) — 0. Foi0 nlLi Hi the singularities of the 
associated ODEs can be written as the solutions of: 

D(4>i, fa) = 0, 

ai • sin(^i) — a n ■ sin(0„) = 0, i = 1, • • • , n — 1, with ^ 
a, -sin(^) = 0, t = 1, N (D.4) 

which are also Landau conditions restricted to pinch singularities of the singularity 
manifolds D(<fii, ipi) = 0. 

Appendix E. The singularities from Landau conditions 

In this Appendix, we give further details corresponding to l|28p . ((29j) obtained from 
the Landau conditions: 

1 - 2w ■ (cos(0j) + cos(^i)) = 0, j=l,---,n, (E.l) 

a.j ■ sin(0j) - a n ■ sin(0 n ) =0, j = 1, • • • , n - 1, (E.2) 

ctj ■ sm(ipj) — a n ■ sin(^„) = 0, j = 1, ■ • ■ , n - 1. (E.3) 

tt 3/i or yf integrand are similar as far as location of singularities of the corresponding ODE is 
concerned. 
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and: 

n n 

<j>j = o, 51 ^ = mod - 2 ^ ( E - 4 ) 

We solve these equations for the values (zero or not) of sin(</> n ) and sm(ip n ). For 
sin(</>„) = sir±(V>n) = 0, the case is simple and gives u> = ±1/4. 

Appendix E.l. The case sin ((/>„) 7^ 0, sin(-0„) = 

In this case, there are k angles ipj — n and the remaining ones are ipj = 0. By lfl7|) . 
the integer k should be even, k = 2 p. From (|E.1[) . we obtain and defin^ 

cos(0+) = J- + 1, cos(<n = ^- - 1. (E.5) 

2w 2w 

One obtains 2p angles <j>j — ±<p + and n — 2p angles cpj — ±cp~ . The angles <pj are then 
partitioned in sets of pi angles +cf) + , {2p — p\) angles —<p + , (n — 2p — p2) angles 
and P2 angles — <p~ . By (|E.4[) . one gets {2p— 2pi)-(p + — (n — 2p—2p2)-<p~. Note that 
some manipulations on the indices lead to cos (|2p| • 4> + ) = cos (|n — 2p — 2k\ ■ <p~) and 
thus \2p\ ■ 4> + = ±\n — 2p — 2k\ ■ <p~ , allowing us to write 

Tip (l/2w + 1) = T„_2p-2fc (l/2w-l), (E.6) 
< p < [n/2], < k < [n/2] - p, 

where T n (x) is the Chebyshev polynomial of the first kind. 

One obtains the same results for the case sin(<^„) = and sm(ip n ) 7^ 0. 



sin(^) = sin(0j) • . , " , j = 1, • • • , n. (E.. 



Appendix E.2. The case sin(<^„) 7^ 0, sin('0 Tl ) 7^ 

In this case, by (|E.2[) . l|E.3|) . we have s'm(<pj) 7^ and sin(ipj) 7^ 0. The equations 
CEp), <(E73|) become: 

cos(ipj) = 1 — 2w ■ cos(<fij), j — 1, • ■ • , n, (E.7) 

sin(V> Ti 
sin( r . 

Squaring both sides of both equations and summing, one obtains 

(cos(0 J ) - cos(^„)) • (cos(<^) - cos(^o)) = 0, (E.9) 

where we have defined 

r , % 4w — cos(</>„) _ 

cos <ft>) = 1 t^t- E.10 

1 - 4w cos(</>„) 

The angles 0j are then partitioned into four sets ±</>o and ±<^„. Note that a 
similar condition l|E.9|) occurs for the angles ipj which are partitioned likewise. Writing 
(|E.T|) . (|E.8P for j = and j = n and with the conditions (fl7|) . the equations become 
in terms of Chebyshev polynomial^: 

(4:W Z \ 

T^Tz) = °' (K11) 

I Note that 4> + and (f>~ (which correspond to ipj = n and ipj = respectively) are not on the same 
footing: indeed, the number of 4> + angles must be even, while the number of <p~ angles depends on 
the parity of n. 

ft Note that in equation llE.lljl one must realise that one takes the numerator of these rational 
expressions. 
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Tni ^2w Z ) Tn2 ^2w 1-Awz 

u ^- i{z) - u ^A^-T^h) 



with 



n\ = p, n 2 = n — p — 2k, (E.12) 

< p < n, < k < [(n-p)/2] (E.13) 

At this step, some computational remarks are in order. In the course of deriving 
(|E.lip . some manipulations such as dividing by a term have been done. Rigorously, 
the solutions that come from (|E,11|) have to be checked against this point. We have 
found, that as they are written, the formulas are "safe" from this perspective, except 
of the following. For n = p/2 (fixing k — for convenience), thus for n even, the 
formulas 1|E.11|) give a common curve which reads: 

w - Irh- (E - 14) 

This relation comes from the condition cos(</> ) = cos(<£„) in (|E.10[) which makes (|E.9|) 
a perfect square. We have checked that considering this condition at the outset, i.e. 
l|E.7j) . (|E.81> yields no solution. 



Appendix E.3. Landau singularities 

We can write the singularities obtained from 1|E.6|) as: 

ra = 3, (l-4w)(l-w)(l + 3w + 4w 2 ) = 0, 

n = 4, (l - I6w 2 ) (1 -W) = 0, 

71 = 5, (1 - Aw) (1 - to) (l + 3 to + 4w 2 ) (1 - 3 to + w 2 ) 

x (l - 7w + 5w 2 -4w 3 ) (l + 8to + 20to 2 + 15t« 3 + 4tt; 4 ) = 0, 
71 = 6, (1 - 16w 2 ) (1 - 4to 2 ) (1 - to 2 ) (1 - 25 w 2 ) (l - 9to 2 ) 

x (l + 3to + Aid 2 ) (l - 3 to + 4to 2 ) = 0. 

The solutions of (|E.llj) include some of the solutions of (|E-6[) . We give in the 
following only those not occurring in (|E.6|) : 

7i = 3, to- (l + 4w)(l + 2to) = 0, 
71 = 4, to = 0, 

7i = 5, to • (l + 4to) (1 + to) (1 + 2 to) (1 + 2 to — 4 to ) 

x (l + 4to + 8to 2 ) (l-to-3to 2 + 4to 3 ) = 0, 
71 = 6, to- (l-10w 2 + 29w 4 ) = 0. 

All these singularities can be identified with the singularities occurring in the 
linear ODE for (n = 3, • • • , 6). 

For n = 7 and n = 8, the solutions of (|E.6[i and IjE.lip can be identified with the 
singularities given in | Appendix C| and obtained in floating point form. They also give: 

ti = 7, (1 + 3to + 4to 2 ) (1 + 4to + 8to 2 ) (1 - to) (l + 2w) 
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x (l-w-3w 2 +4w 3 ) = 0, 
n = 8, (1 -3u) + 4iu 2 ) (1 + 3w + Aw 2 ) (l - 10 w 2 + 29w 4 ) = 0, 
which have not been found in the series with the currently available number of terms. 

Appendix F. Heegner numbers and other selected values of the modular 
j-function 

The nine Heegner numbers [51] and their associated modular j-function j(r), yield 
the following conditions in the variable w: 

+ = (12) 3 , (l-8w 2 ) (l - 16u> 2 -8u> 4 ) = 0, 

j(l + iV2) = (20) 3 , (64 w 4 + 16 w 2 - l) 

x (64 w 8 + 1792w 6 - 368w 4 +32w 2 - l) = 0, 

jfl+^j = (0)3, l-16u- 2 + 16 W 4 = 0, 

j(i±M) = (_i5)3 j (1 -31 W 2 + 256^ 4 ) (1 ~16^ 2 + W 4 ) 

x (1 +3w +4w 2 )(l -3w + 4w 2 ) = 0, 



and 



. /I + iVd\ 
\—2-J 



= (-32) 3 , P 3 = l-48w 2 + 816w 4 - 5632w 6 
+ 45824 m; 8 - 536576 w 10 + 4096 w 12 = 0, 



(-to) 3 , P d = with: 

P d = P 3 + N- (l - 16 w 2 ) ■ w 8 , 
with the following values for the triplet (d, m, N): 

(19, 96, 851968), (43, 960, 884703232), 

(67, 5280, 147197919232), (163, 640320, 262537412640735232) 



Beyond Heegner numbers there are many other selected quadratic values [52, 53J 
of j, for instance: 

3 = -4096- (l5 + 7^ 3 = jC- 1 ^) (F.l) 



2 

Which is known [51] to be one of the eighteen numbers having class number h{—d) = 2, 
and which corresponds to the quadratic relation -134217728000 + 117964800 j+j 2 = 
0. Recalling the expression of the modular j-function in term of the variable w 

(1 - 16w 2 + 16w 4 ) 3 
j = (1 - 16 w 2 ) w s ' (F ' 2) 

this quadratic relation in j becomes a quite involved polynomial expression that we 
have not seen emerging as singularities of (the linear ODE's of) our n-fold integrals. 
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Appendix G. Landen transformations and the modular j-function 

In this Appendix the modular j-function (|53| will be seen, alternatively, as a function 
of the modulus k, and thus denoted j[k], or as a function of the half period ratio r, 
and thus denoted j(r). The modular function called the j-function when seen as a 
function of the modulus k reads: 

(l-fc 2 + £: 4 ) 3 

(G ' 1) 

Increasing the modulus by i|48p . the modular function j(k) becomes: 

JM-J.W- 16. ( 1 + 14fc3 + fc ;)' . (G.2> 

Iterating this procedure once more one obtains: 

(k 4 + 60k 3 +134/c 2 + 60fc + l) 3 

^ = *>w = 4 ■ - — fc . ( i +fc ) 2 (i-fe) s ■ (G - 3) 

The decrease of the modulus by fl47| gives: 

r , r , (fc 4 -16fc 2 + 16) 3 , 

j[k-i] = j-i[k] = 16 " fc8(1 _ fc2) J ■ (G- 4 ) 

The next iterations (the cube of (|i"8jl and the square of lj47| ) gives algebraic 
expressions for j[k]. 

It is easy to get a representation of the Landen transformation on the modular 
j -functions by elimination of the modulus k between lj53"l) and (|G.2p . One obtains the 
well-known fundamental modular curve |4~9l 1501: 



I'M./- h) = f ■ 3i - U +h)- {j 2 + 1487 331 + 3 2 i) 

+ 3 ■ 15 3 • (16 j 2 - 4027 3 3i + 16 j?) (G.5) 
- 12-30 6 -C/ + 8-30 9 = 0. 

This algebraic curve is symmetric in j and ji . We will obtain the same modular 
curve 1|G,5|1 by elimination of the modulus k between (|G.2j) and (|G.3[) . or between 
l|G.l|l and l|G,4|l . The two modular functions j and j± are invariant by the SL(2, Z) 
modular group, and, in particular, transformation r — > 1/r. As a consequence, the 
transformation r — > 2-r, and its inverse r — > r/2, /ia«e to be on the same footing in 
the modular curve representation (|G.5j) for the Landen and Gauss transformations. 

Similarly, one can easily find the (genus zero) modular curve T2 obtained by the 
elimination of the modulus k between (|G.ip and 1|G,3|1 . (or between (jG.4|) and (|G.4P ). 
which corresponds to the transformation r — * 4 • t and, i/*e same iime, to its 
inverse t — > r/4. This last algebraic curve is, of course, a/so a modular curve. 



Appendix G.l. Fixed points of these modular representations in terms of j-function 

Transformations like j —* ji, or j — > 32, corresponding to the previous modular 
curves, are not (one-to-one) mappings, they are called "correspondence" by Veselov O 
[24] . In order to look at the fixed points of the Landen, Gauss transformations (or 
their iterates) seen as transformations on complex variables, within the framework 
of (modular) representations on the modular j-functions, we write, respectively, 
IiO'jji = 3) = and T 2 (j,32 =j)=0 
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The "fixed points" Ti(j,ji = j) = of the (modular) "correspondence" (|G.5[) . 
are j = j x = (12) 3 or (20) 3 or (-15) 3 . 

The "fixed points" T 2 {j, j% = j) = of modular curve corresponding to the square 
of the Landen transformation, are j = — (66) 3 , or 2 • (30) 3 , or (— 15) 3 or the 
solution^ of j 2 + 191025 • j - 121287375 = 0, namely: 

• (5 + 4 • v 5) = j[r = (G.6) 



3 ~ V 2 
and its Galois conjugate (change \fh into — y5) 



Appendix G.2. Alternative approach to fixed points of the Landen transformation and 
its iterates 

In order to get the "fixed points" of the Landen transformation, let us impose that 
(|G.1|) and (|G.2j) are actually equal, thus j[k] = j[k{\. This yields the condition 
(already seen to correspond to the x -singularities 1 + 3w + 4w 2 = 0): 

(4 k 2 + 3 k + 1) (fc 2 + 3 k + 4) = (G.7) 

together with: 

(4fc 2 -3fc+l) (fc 2 -3fc + 4) 

(fc 2 + 2 fc - 1) (fc 2 - 2 fc - 1) (1 + fc 2 ) = 0. (G.8) 

The first two polynomial conditions in 1|G.8|) . (4 fc 2 - 3 fc + 1) (fc 2 - 3 k + 4) = 0, 
correspond to the Heegner number associated with the integer value j = (— 15) 3 . 
The next two polynomial conditions in 1|G,8|1 . fc 2 ±2fc— 1 = 0, correspond to the 
Heegner number associated with the integer value j = (20) 3 . The last polynomial 
condition in 1|G.8|) . 1 + fc 2 = 0, corresponds to the Heegner number associated with 
the integer value j — (12) 3 . 

Similarly, in order to get the "fixed points" of the square of the Landen 
transformation, let us require that (|G.1[) and (|G.3[) are actually equal: j[k] — jfa]- 
This yields the conditions (|G.7j) (fixed points of the Landen transformation) together 
with: 

(fc 2 -6fc+l) (l + 14fc 2 + fc 4 ) = (G.9) 

(fc 4 -6fc 3 + 17fc 2 + 36fc + 16) (G.10) 

x (16 fc 4 + 36 fc 3 + 17 fc 2 - 6 k + 1) = 0. 

In (|G.9|1 the condition 1 + 14 fc 2 + fc 4 = (or 1 — 16w 2 + 256 w 4 = 0) corresponds 
to j = 2 (30) 3 which is not a Heegner number but actually corresponds to complex 
multiplication. The condition fc 2 — 6fc+l = in 1|G.9|) (or 1— 32 w 2 = 0) corresponds 
to j = (66) 3 which is not a Heegner number either but actually corresponds to 
complex multiplication. Note that both polynomials under the Landen transformation 
|48l) give respectively j = (0) 3 and j = (12) 3 , i.e. Heegner numbers. The last two 
(self-dual) conditions in (|G.10|) . read in w 

1 -9w +I7w 2 +24w 3 + 6w 4 = 0, (G.ll) 
1 + 9 w + 17 w 2 - 24 w 3 + 6 w 4 = 

and yield as selected value [52, 53] of j, the quadratic roots —121287375 + 191025 j + 
j 2 = 0, already given in 1|G,6|) . 



| This corresponds to a value of j of class number h(—d) = 2, see (58) in |51| . 
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One more step can be performed writing the condition i[fc-i] = One gets 

the conditions: 

(fc 2 + 3 fc + 4) 2 (4 fc 2 - 3 fc + 1) (fc 2 + 2 fc - l) (fc 2 + l) = 

previously obtained and corresponding to j = (— 15) 3 , 20 3 , 12 3 , together with: 

fc 6 - 27 /fc 5 + 363 fc 4 + 423 fc 3 - 168 fc 2 - 144 fc + 64 = 0, 

fc 6 + 17fc 5 + 143fc 4 + 203fc 3 + 52fc 2 + 32fc + 64 = (G.12) 

corresponding, respectively, to the two cubic relations on j: 

1566028350940383 - 58682638134 j + 39491307 j 2 + j 3 = 0, 

12771880859375 - 5151296875 j + 3491750 j 2 + j 3 = 0. (G.13) 

These conditions (|G.13|1 yield quite involved polynomial expressions in the variable w 
that we have not seen emerging as singularities of (the linear ODE's of) our n-fold 
integrals (or the Y^ n ' or $( n > either). 

Appendix H. Linear differential operators for the Sorokin integrals 

Recall the occurrence of zeta functions evaluated at integer values in many n-fold 
integrals corresponding to particle physics, field theory, ... For instance, the following 
integral [33, 36J is associated with C(3): 

f 1 (l-u) n u n -(l-v) n v n -(l-w) n w n 
Ijz) = / dudvdw — (H.l 

J [1 — u v) n+L ■ (z — uv w) n+L 

From the series expansion of this holonomic n-fold integral, we have obtained 
the corresponding order four Fuchsian linear differential equation. On these linear 
differential operators the "logarithmic" nature of these integrals becomes clear. 

The fully integrated series expansion of the triple integral IjH.ip is given by (where 
x denotes 1/z): 



DC 



n+i+ 



1 r 2 (n + i)-r 4 (n + i + i) 

T(i+ 1) -r 3 (2 + 2n + i) 



3 F 2 {n + 1, n + i + 1, n + i + 1; In + i + 2, 2n + % + 2; 1). 

The triple integral I n (x) is solution of the order four Fuchsian linear differential 
operator (Dx denotes d/dx) 



■Dx 2 



[x — 1) x 

(7 x 2 + (n 2 + n-5)x - 2n(n + l)) 

(x — l) 2 x 2 

(x 2 + 2 n in + 1)) 

1 l LL . d x 

(x — l) 2 x 3 

n (n + 1) • ((n 2 + n + 1) x + (n - 1) (n + 2)) 



(x- l) 2 x 4 

which has the following factorization 



Dx + d -^l) ■ (Dx + ^1) (H.2) 

dx J \ dx J 

x (Dx + d -^) ■ (Dx ■ dHM) 



dx ) V dx 
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where: 

A x = -(n - 1) • ln(ar) + 2 • ln(x - 1) + ln(P„), 
A 2 = (n + 1) ■ In(ar) - (n - 1) ■ In(a; - 1) - ln(P„) + ln(Q n ), 
A3 = -n ■ ln(x) + (n + 1) ■ ln(x - 1) + ln(P„) - ln(Q„), 
A 4 = ra-ln(x) -ln(P„), 

and where P n and Q„ are polynomials in x of degree n. They are the polynomial 
solutions behaving as • • • + x n for a system of coupled differential equations (pl™^ 
(resp. Qn™"*) denotes the m-th derivative of P„(x) (resp. P„(x)) with respect to x): 

(x - I) 2 • x 2 • PW 

- 2 • (2 (x - 1) • n - 3x + 1) • (x - 1) • x • P r { 3) 

+ ((2x-l)(3x-4)-n 2 - (l2x 2 -13x + 2) - n +(7x-5) x) • P^ 2) 

- (2 (2a; -3) • n 3 - 2 (3 x - 1) • n 2 + 2 (2 x - 1) • n - x) • P,W 
+ n 4 -P„ = 0, 

- (x - 1) • x • P„ • 

+ (2 (x - 1) • x • P« + (1 - x + 2n) • P„) • Q« 

- (2 (x - 1) • x • P^) - 2 ((x - 2) n - x) ■ P^ + n 2 P„) ■ Q n = 0. 
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